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Density profiles of small Dirac operator eigenvalues for two color QCD at 
nonzero chemical potential compared to matrix models 
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We investigate the eigenvalue spectrum of the staggered Dirac matrix in two color QCD at finite chemical 
potential. The profiles of complex eigenvalues close to the origin are compared to a complex generalization of the 
chiral Gaussian Symplectic Ensemble, confirming its predictions for weak and strong non-Hermiticity. They differ 
from the QCD symmetry class with three colors by a level repulsion from both the real and imaginary axis. 



1. Introduction 

The comparison of small Dirac operator eigen- 
values from lattice QCD to matrix models has 
been very successful at zero chemical potential 
At fj, ^ detailed comparisons to the microscopic 
Dirac operator origin spectrum have been few 
for two reasons. First, unquenched QCD suffers 
from a complex action problem and second, ma- 
trix model predictions where not available. Very 
recently substantial progress has been made for 
complex chiral matrix models. Quenched |2l3j 
and unquenched 0] predictions for QCD have 
been achieved, as well as results for adjoint 
gauge theories and staggered two color QCD 5 . 
They have been successfully tested for quenched 
QCD and our aim is to extend these results to 
other symmetry classes. The advantage of SU(2) 
is to avoid the complex determinant of the action 
matrix |7]. Already a few years ago, we computed 
the lowest eigenvalues of the unquenched SU(2) 
Dirac operator at 7^ jH], with low statistics 
though in). We redo the analysis as we can now 
compare to jH]. For simplicity we only present 
quenched results here. 



2. Matrix model predictions 

In Ref. |10| a complex chiral matrix model with 
the same global symmetries as QCD at /i 7^ was 
introduced. Its microscopic Dirac spectrum was 
obtained only very recently j3l4j . after previous 
asymptotic results in [2] for a related model. Al- 
though the presence of fi alters the chiral symme- 
try breaking patterns for SU(2) and the adjoint 
representation, the symmetry analysis of Dirac 
matrix elements being real, complex or quater- 
nion remains valid ,41, . Following the two-matrix 
model approach ^ new results for quaternionic 
matrices are now available |S] . We only reproduce 
the quenched microscopic density here and refer 
to Ref. |Sj for more details. 

In QCD at /X = 0, the lowest Dirac eigenval- 
ues scale with the volume 1/1^ to build up a fi- 
nite condensate S according to the Banks-Casher 
relation. At weak non-Hermiticity, a concept 
first introduced in jl2) . this scaling remains un- 
changed [^. The support of the density remains 
quasi one-dimensional as we send fi ^ 0, keeping 
limAT^oo 2iV^^ = fixed. Here 2^ is the 

size of the random matrix. In contrast to that, at 
strong non-Hermiticity the eigenvalues fill a two- 
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dimensional surface, and thus the scaling is mod- 
ified to ^ 1/VV 2 . The quenched microscopic 
density at weak non-Hermiticity in the sector of 
topological charge i' is reading [31 

X ds dt t-5e-^(i+*)"' 
Jo Jo 

X (j2.(2V^oj2.(2^/ir)-(e^r))- 

Here, we have rescaled the complex eigenvalues 
according to ^ Vz. 

The limit of strong non-Hermiticity can be ob- 
tained by taking the limit a ^ oo in Eq. The 
resulting microscopic density is symmetric with 
respect to the axes 3fie(^) and 9m(^). 

The zero momentum sector of chiral pertur- 
bation theory with /i ^ contains two param- 
eters ^21) ^ and the decay constant with the 
relative rescaling between mass and jj. term de- 
pending on the ratio only. In the weak limit Eq. 

becomes a function of ^/a only, by substitut- 
ing sa^ = r. Note that the first integral then 
runs to . However, at small a > this change 
of variables breaks down. The integral becomes 
approximately a-independent while the prefactor 
strongly varies with a. Thus in this regime eigen- 
values z (or masses) and jjL can be rescaled inde- 
pendently with two parameters. 

3. Lattice data at /i 7^ 

Our computations with gauge group SU(2) on 
a 6-* lattice at /? = 4/g2 = 1.3 have Nf ^ 2 fla- 
vors of staggered fermions. For this system the 
fermion determinant is real and lattice simula- 
tions of the full theory with chemical potential are 
feasible. To compare with the quenched predic- 
tions from above we have chosen a large value for 
the quark mass ma = 20 (a being the lattice spac- 
ing in physical units) that effectively quenches the 
theory. We have checked this by comparing to 
simulations with ma = 2, following exactly the 
same analytical curve as in Fig.nbelow. A phase 
transition occurs at /ic « m^/2 « 0.3 where the 
chiral condensate goes to zero and a diquark con- 
densate develops [7]. 



At strong gauge coupling and thus away from 
the continuum limit, staggered fermions have the 
disadvantage of shifting the (topological) Dirac 
zero modes and mixing them with the nonzero 
modes. We have accounted for this by setting 
z/ = in Eq. above. Only more sophisti- 
cated methods as improved actions [2] can re- 
cover topology in the staggered approach. We 
produced around 10000 configurations for two 
values oi n — 0.001 and ^ — 0.2, correspond- 
ing to weakly and strongly non-Hermitian lattice 
Dirac operators, respectively. 

Fig. n shows that the data agree very well 
to the functional form of Eq. |^ at weak non- 
Hermiticity. After normalizing histograms and 
formula to unity we rescale the lattice eigenval- 
ues z hy ^ = nz/d, with the spacing d oc 1/V. 
The weak non-Hermiticity parameter has been 
rescaled with a second parameter ^ f^/d to 
match the data. In contrast to QCD with three 
colors |2I3I4| the symplectic symmetry class shows 
an additional suppression of eigenvalues along the 
real axis, resulting into the double peak struc- 
ture. The zero modes are driven away from the 
real axis for /i as small as /i = 0.000001. This 
indicates that the limit /i ^ is discontinuous. 

Fig. |2] shows data for ma — 20 at strong non- 
Hermiticity, with ^ = 0.2 still below the phase 
transition. Here, the eigenvalues are rescaled with 
the square root of the volume, ^ '-^ z/Vd. Lacking 
an analytical prediction for the level spacing we 
have rescaled our eigenvalues by one parameter to 
fit Eq. The normalized data nicely confirm 
the prediction, showing the correct symmetry. 

4. Conclusions 

We have tested the analytical predictions from 
a complex chiral symplectic matrix model using 
effectively quenched two color QCD. The agree- 
ment demonstrates that also at fj, ^ matrix 
models are an excellent tool to describe the dif- 
ferent gauge theory symmetry classes. In this 
preliminary study we have established the cor- 
rect qualitative behavior from quenched data. A 
detailed analysis including different lattice vol- 
umes, scaling, and dynamical flavors with small 
mass is in progress. 
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Figure 1. Data (histograms) for V — 6"^ at fi = 
0.001 vs. Eq. for weak non-Hermiticity (a ~ 
0.4), cut parallel to the real axis along maxima 
(upper plot) and to the imaginary axis at the first 
maximum (lower plot). 
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